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The (discrete) Gross-Neveu model is studied in a lattice realization with an ,/V-component Majo- 
rana Wilson fermion field. It has an internal O(A0 symmetry in addition to the euclidean lattice 
symmetries. The discrete chiral symmetry for vanishing mass is expected to emerge in the con- 
tinuum limit only. The lattice theory is first recast in terms of two-valued bosonic link variables 
(dimers). In this representation, which coincides with the loop representation obtained earlier 
by Gattringer with the help of eight-vertex-models, the Boltzmann weight is essentially positive. 
While standard local updates are possible in this form we construct a further exact transformation 
where we generate dimer configurations as Peierls contours of an Ising model with a local action 
residing on plaquettes. For this model a Swendsen-Wang type cluster algorithm is constructed. At 
vanishing coupling it is numerically demonstrated to almost completely eliminate critical slowing 
down. Although further tests are required, an avenue to numerical studies of the Gross-Neveu 
model with unprecedented precision seems open. 
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1. Introduction 

The O(A0 symmetric Gross-Neveu (GN) model [JTJ] forms a particularly simple family of self- 
interacting fermionic quantum field theories. It may be regarded as the natural fermionic counter- 
part of the [Af-component] 4 theory of real bosons. While they are renormalizable in four dimen- 
sions the same is true in d = 2 for GN, where it is even asymptotically free. It is most naturally 
written in terms of Majorana fermions t, which are pairwise related to the more familiar Dirac 
fermions by 

y= -^(fr + ife), \jf=-j=(tf-i£)tf, ^ = -^ T : charge conjugation. (1.1) 
We here define the action immediately on the lattice with Wilson-fermions as 

We work with lattice units, a = 1. The first term contains the mass and the diagonal part of the 
Wilson term and is followed by the interaction and the hopping term. Note that for Majorana 
fermions the backward hopping terms coincide the forward ones. All i^-bilinears are contracted 
in the internal N-valued index making the 0(/V) symmetry manifest. It allows only one type of 
4-fermi interaction, hence the strict renormalizability. The discrete chiral transformation £, — > 75^ 
is broken only by the mass and the Wilson term. As with chiral symmetry in QCD it is expected 
to emerge in the massless continuum limit where the bare Wilson mass has to be tuned to a critical 
value m c . 

The partition function follows, 

,2 7)2 



/ D ^ = exp{^£^}(3,[».M])~. (1-3) 



By allowing for an ^-dependent mass m we re-wrote it here in terms of the partition function Z ( 







of one Majorana fermion in such a background. Its action is the bilinear part of (1.2). In ref.[g] 
more details on everything said here can be found. In particular it is explained, how the interacting 
theory can be simulated once one has the efficient algorithm for Zo[wi(jc)] available that we study 
below. Such simulations are presently under way. 

2. Transformation to a dimer ensemble 

To expand the partition function with <p(x) = 2 + m(x) we first exploit the nilpotency of even 
Grassmann elements and obtain. 

Z = jD^Ylil + cp^Yl^ + ^ix^-^^x + fi)^. (2.1) 

Next the product of the hopping terms on each link is organized with the help of link dimer- 
variables k(x,jj.) = 0, 1, 

z o= 1 /^n^+^^nf^w^^s^+A))* 1 * ( 2 - 2 ) 

{k{x,n)} J x x,H \ Z J 

= £ p[*]xsign[*] (2.3) 
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Figure 1: Lines correspond to dimer loops which also form contours separating domains of up and down 
spins. The four different possible loop topologies on the torus are exemplified. 



with 

<p(x) ifnodimeratx(= monomer) 

t-t I 1 if2dimers in the same direction at x _ , 

p[k] = TT x < . r . (2.4) 

11 1 l/v2 if 2 dimers in different directions at x 

else 

To arrive at the weight p one has to observe that the Grassmann integral constrains dimers (links 
with k(x,n) = 1) to form closed loops and at each site either two or zero dimers can touch. A site 
with no dimer adjacent is called a monomer and leads to a factor (p(x). The weight for such loops 



can be computed analytically [|5|] and leads to ( |2.4[ ). The factor sign[fc] equals +1 unless dimer 



loops wind around the boundary, in which case it depends on whether the boundary conditions are 
periodic or antiperiodic. This leads to relations between fermion and dimer-ensembles of the type 



=Zf ) +Z l k " -Z^+Zl 1 (+3 more such relations). (2.5) 
In this formula is a fermion (anti)periodic in (time) space, and Z\ Q is a partition function of the 



type (2.3) restricted to configurations with 1 (mod 2) loops around the torus in time and in space 
and analogously for the other cases. The four possibilities are illustrated by the lines in Fig.[jJ An 
example of the translation of observables between the representations, derived by taking derivatives 
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Figure 2: Autocorrelation times for local (x) and cluster (*) algorithm for a massless fermion. 



with respect to m(x), relates the scalar fermion density and the monomer density 
2 + m 



I fermion 



1 if monomer at x 
else 



(2.6) 



3. Spin representation and cluster algorithm 

In the form of the dimer ensembles we have recast the originally fermionic theory in terms of 
discrete commuting variables with local interactions. They can be updated by a local Metropolis 
algorithm as demonstrated in [ft]. Due to the constraints [the zero in ( |2.4] )] the minimal updates 
have to change dimers around a plaquette. Such updates do not change the dimer loop topology 
which is determined by the initial state. Such a fixed topology dimer ensemble corresponds to a 
combination of fermion ensembles with several kinds of boundary conditions. The local simulation 
of such an ensemble has roughly the same complexity as simulating a standard Ising model. 

Also as in the Ising model the local update leads to critical slowing down with a dynamical 
exponent z around two. For the case equivalent to a free Wilson Majorana fermion at vanishing 
mass, which is critical, the upper crosses in Fig|2]at least qualitatively show the expected growth of 
the integrated autocorrelation time for the monomer density. In the Ising model this severe deterio- 
ration of efficiency is practically completely eliminated by collective cluster update procedures like 
the one of Swendsen and Wang []6p. The good news is now that also here, for our Ising type theory 
with plaquette interaction - exactly equivalent to Wilson fermions - we have found an equally ef- 
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ficient technique. We briefly sketch it here, but for a more detailed understanding probably [Q] has 
to be consulted. 

We consider the lattice dual to the one carrying the dimers (and originally the fermions). Its 
sites can be drawn in the centers of the original plaquettes. This is where we put Ising spins, which 
appear as plus and minus signs in Fig.|l]. We also read off there, how the closed (non-intersecting) 
dimer loops are obtained as the Peierls contours encircling areas of a given spin orientation. As dis- 
cussed in [^] this is an exact two-to-one mapping of Ising configurations obeying a certain plaquette 
constraint and allowed dimer configurations. The four different dimer topologies correspond to the 
four possible combinations of (anti)periodic boundary conditions of the Ising spins. There an- 
tiperiodicity forces domain boundaries into our 'magnet' which are precisely the non-contractable 
dimer loops, see again Fig.[l|. The Ising partition function - for one choice of topology/boundary 
conditions — is now given by 

izf=zt-= i n AtX (3 - 1} 

{s(x)} plaquettes \ A 1 ^ ) 

The desired weights for the individual configurations dictate the values of w, 

w(l l)=<p(x), W (lr)=-L, w ( — \ =lf w (- ) (3.2) 



with other values obtained form these by rotating or flipping the four spins. The essential observa- 
tion for the cluster approach is, that this function of four spins around a plaquette may be written 
as a linear combination of ten terms consisting of Kronecker deltas like S S2 iS3 tying together pairs 
of spins with positive coefficients, 



w 



+q( H I +r U + l H KH- (3-3) 



In the above symbolic formula the delta-bonds are represented by the solid lines. By considering 
all 16 configurations of the four spins (and the obvious symmetries) the coefficients 

r=^M p = J—- r tf = i_J_ + r (all ^ Oif m < 2\fl) (3.4) 

result. This is a rather direct generalization of the Fortuin-Kasteleyn bonds in where only just 
pairs of spins have to be considered. We now introduce correspondingly 10-valued variables b 
on the sites of the original (=plaquettes of the dual) lattice as additional variables. The partition 
function is now given by a sum over both the spin and the new bond field, 



Z s = I n^M-vAM • (3-5) 

{2>(.v) = l,...,10;s( 2 )=±l}plaq \ J 1 S 2 ) 

Here the factors A S 0, 1 depend on both bonds and spins and contain the 8 associated with a given 



b- value, while the P are the weights from (3.4). The decisive trick is now to refrain from exactly 
summing over b, but to alternatingly do Monte Carlo updates to both kinds of variables. In a first 
pass, for fixed spins, at each site a new bond is picked among those allowed by A with relative 
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probabilities given by P. This amounts to a local heatbath, and the ^-updates at different sites are 
independent in this step. With new bonds given we next pick a new spin configuration. All spin 
configurations on the whole lattice in this step are either forbidden by A = or have the same 
positive weight. If one chooses one of the latter at random with equal probability this furnishes a 
global heatbath in this step. With the help of percolation cluster search algorithms which divide the 
spins into groups ('clusters') that are tied together by the 8 bonds, this becomes possible at a cost 
of O(volume) only. This is the nonlocal update step that (almost) eliminates critical slowing down 
as witnessed by the stars in Fig.||. 

In [^] two important generalizations are given that we only mention here. First, in the step 
of picking new spins at fixed bonds one may also at the same time consider to switch among the 
four possible boundary conditions. One then simulates with 'fluctuating' boundary conditions the 
ensemble Z+ + +Z~ + +Z+~ +Z S . By reweighting with a factor depending on the boundary con- 
dition as a dynamical variable one can then also effectively measure at fixed fermion boundary con- 
ditions, see ( |2.5| ). This factor is of fluctuating sign which is a remnant of the fermionic minus sign 
problem, still visible at finite volume here. It leads however to no severe cancellations that cannot 
be handled. Moreover it even seems possible to construct an improved estimator at least for some 
quantities in the reweighted ensemble. The second generalization allows for (p(x) = 2 + m(x) < 2, 
i.e. negative mass. If the simulations studied for a fermion in the 'external' field m(x) here are ac- 
tually embedded to simulate the interacting GN model, such local masses cannot be avoided. The 
decomposition (33) with weights (3_4) requires m>0. There is however another decomposition 
available in the case of negative mass. It contains both bonds as before and antibonds like 8 Sl - S2 
with then altogether 14 terms and positive weights again. Clusters now also contain tied-together 
opposite spins, but otherwise the construction of clusters which are flipped collectively with prob- 
ability 1/2 is very similar. With local m(x) some plaquettes may use the one and some the other 
decomposition. 

Both in [§] and here only free fermions have been simulated numerically. In the dimer or Ising 
form this case in fact does not seem so very special, but of course from the fermionic realization 
exact results are available for all quantities. In the meantime first simulations are running however 
with interacting N = 2 and TV = 8 fermion flavors. First results to be compared with [Q] look con- 
sistent and the update seems to remain very efficient. Unfortunately the method so far is restricted 
to two dimensions. Fermions, and the possibility to 'bosonize' them, is known to be special to this 
dimensionality. Nevertheless, it seems worth to keep thinking about 'non-determinant' approaches 
to fermions also for d = 3,4. 

Acknowledgement: The author thanks the Deutsche Forschungsgemeinschaft (DFG) for support in 
the framework of SFB Transregio 9. 
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